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It is shown that the eigenstates of the quantized simplest closed (elliptic) Nambu 
- Goto string, called glueballs, have quantum numbers I G j PC = + j ++ . Lightest 
glueballs have spins j = 0, 1 and 2 and the same mass 1500 ± 20 MeV. They 
correspond to /o(1500), /i(1510) and /2(1565)- mesons. Next glueballs have j = 
0,1,2,3,4 and the same mass 2610 ± 20 MeV. The slope of the glueball Regge 
trajectories is twice as small as for qq-mesons. The intersept of the leading glueball 
trajectory — the pomeron Regge trajectory — is 1.07 ± 0.03. 



1 Introduction 

Overwhelming majority of known mesons can be described as eigenstates of 
the quantized straight-line Nambu - Goto string with the Dirac quarks at 
the string ends [1]. The string approximately accounts for the nonpertur- 
bative contribution of the gluon field at large distances and is characterized 
by a string tension parameter. The quarks are characterized by their cur- 
rent masses. All theese parameters, together with phenomenological param- 
eters describing nonstring short-range gluon contribution, are determined 
from comparison of the theoretical and experimenral meson mass spectra for 
mesons lying on the leading Regge trajectories. 

Since this approach has appeared to be successful for the quark-antiquark 
mesons, it seems natural to generalized it to the simplest closed strings with 
the aim to describe the gluon meson states called glueballs. 

The classical and quantum description of the simplest closed (elliptic) 
Nambu - Goto string was done in Ref . [2] . In the present paper a nonstring 
short-range gluon contribution is introduced and possible iterpretation of the 
results of Ref. [2] is considerd together with the qq-mesoia analysis of Ref. [1]. 

In Sec. 2 the results of Ref. [2] are obtained by the 1-form method of Ref. 
[3] which is simpler than the Dirac brackets method of Ref. [2]. Properties 
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of the classical solution, its stability in particular, are also considered. The 
solution is not stable when the string mass is zero. This value of mass must 
be excluded when quantizing the string (otherwise the quantization would 
lead to a tachion). In Sec. 3 the canonical quantization, glueball wave 
functions and glueball Regge trajectories are considered. Comparison with 
experiment allows one to identify the lowest glueball states, to fix the only 
free parameter of the model and to calculate the model predictions: masses 
and quantum numbers of higher glueballs and the pomeron Regge trajectory. 
In conclusion, the obtained results are briefly discussed. 

2 Classical glueball model 

A closed string is described by a 4- vector r), which depends on relativis- 
tic scalars characterizing position of points on the string (er) and evolution 
of the string (r); < a < 2rc, rr M (0,r) = x M (27r,r). Prime and dot denote 
partial derivatives with respect to er and r, respectively. The Nambu - Goto 
Lagrangian 

1-2-K 

C = -a ((x'xf -x n x 2 ) 1 ' 2 da, (1) 
Jo 

where the string tension parameter a is known from the qq- meson analysis 
[1], can be represented in the Hamiltonian form [4] 

p2tt p2tt 

£ = - pxda- (hAa 2 x' 2 + p 2 ) + h 2 px') da. (2) 
Jo Jo 

Here the first term determines the Poisson brackets of the components of the 
string coordinate x and the conjugate momentum density p [3] and the second 
term is the string Hamiltonian which, due to the Lagrangian symmetry, is a 
linear combination of the constraint functions (hi are Lagrange multipliers). 
Let us consider the simplest (elliptic) cofiguration of the closed string 

x(a, r) = r(r) + ?i(r) coscr + (72 (r) sino" (3) 

and use parametrization in which the momentum density has the same con- 
figuration 

p(a, t) = 7r _1 (2 _1 P(r) + 7Ti(t) cose + 7t 2 (t) sine). (4) 
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The orthonormal parametrization, for instance, in which p = ax, or parameteriza- 
tion in which x'x = and x 2 /x' 2 does not depend on a, possess this property. 

Putting (3) and (4) into (2), we get constraints which correspond both to 
the Lagrangian symmetry and to the choice of parametrization 

10 

C = -Pr - 7ri<ji - vr 2 g 2 - c ^> ( 5 ) 

i=i 

where q are proportional to integrals over da of /ii j2 , multiplied by 1, cos ka 
or sin ka, k = 1 or 2. The constraints 0j are 

1 = aV(g 1 2 + g 2 2 )+2- 1 P 2 + 7r 1 2 + 7 r 2 2 ) (6) 

02 = VTig 2 - 7T2?1 (7) 
03 = Pm, 04 = P?f2, 05 = Pq2, 06 = ^<7l, (8) 

07 = a 2 vr 2 (g 2 - q 2 ) + tt 2 - vr 2 , 8 = -a 2 ix 2 q x q 2 + 7Ti7r 2 , (9) 

09 = VTig 2 + 7T 2 gi, 010 = 7T2g2 - 7Tl9l- (10) 

The 1-form in (5) determines the Poisson brackets of the introduced variables. 
The nonzero brackets are 

{pv"} = K,tf} = *r, < = i,2. (ii) 

The constraints (8), (9), (10) are second kind constraints with respect to 
theese brackets (due to the choice of parametrization) and must be solved ex- 
plicitly. To this end we first consider conserved and parametrization-invariant 
string variables which have zero brackets with the Hamiltonian in (5). Such 
variables are the total string momentum P and its nass m = \f~P 2 and the 
string spin 

J, = E W^V<, v v = P v /m. (12) 

i=l,2 

J M has zero brackets with P v and with string Lorentz scalars. The Poisson 
brackets of the spin with string (pseudo) vectors Y are 

a Y^. (13) 

It is remarkable that the elliptic string has one more conserved parametrization- 
invariant pseudovector [2] 

W, = e^v^ian)- 1 ^ + antfqS). (14) 
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Its brackets with all constraints (6 - 10) vanish and 

{W II ,W u } = e IIua pv a jP. (15) 
Let us introduce so-called string pseudospins 

L 12 = 2~ 1 {J±W). (16) 

Their brackets are 

{Lin, Ljv] — Sij^fj,ual3V a L^ , (17) 

{L t)l ,P l/ } = 0. (18) 

Thus, the elliptic string has two independent conserved pseudospins the sum 
of which is the string spin 

J = L 1 + L 2 . (19) 

Now we can write down the solution of the constraints (8), (9) and (10). Let 
us introduce the tetrad of vectors e£(P), a = 0, a; a = 1, 2, 3 

e a ep = g a p, e = v, e^pv e o e a e b e c = e abc (20) 

and expand the vectors Y = q iy 7Tj, L iy J, orthogonal to P, with respect to the 
tetrad 

Y» = e»Y a , a =1,2,3. (21) 
We shall use the notations of 3-vectors for a set {Y a } 

{Y a } = Y, {e abc n b m c }=[n,m} (22) 

and so on. Then the solution of the constraints is 

qi = -^([Li,n 1 ] + [L 2 ,n 2 ]), tti = — (ni + n 2 ), (23) 

<& = t— («i - ^2), vf 2 = ^([Li,ni] - [L 2 ,n 2 ]), (24) 
47ra 4 

where 

= I I, (25) 
ra? = 1, = 0. (26) 
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Instead of 16 variables <zf,7rf with 8 constraints, we have introduced 12 vari- 
ables Lf,nf with 4 constraints (26), which can be easily solved, for instance, 
using spherical coordinate systems. 

Putting the solution (23 - 26) into the Lagrangian (5), we get 

C = -PZ+ J2 v / 5-^)-c2(v / 5-v / 5)> (27) 

i=l,2 47rfl 

where the coordinate conjugate to the total momentum is equal to 

Z» = r» + \e abc e^J c . (28) 
Nonzero brackets following from the 1-form in (27) are 

{P",Z"} = <r, (29) 

{Lia, Lib} — tabcLic, {Lj a , n^} = tabetic- (30) 

The constraints in (27) and the brackets (29), (30) were first obtained in Ref. 
[2] by a different method. 

Using stationary condition for (27), it is not difficult to obtain the classical 
solution for the elliptic string. In the orthonormal parametrization (c\ = 
1, C2 = after variations) 

x(a,r) =x + -^-Pr + e a q a ((r,T), (31) 
27ia 

q(a,r) = - — (— sino;sin((T + a ) cos(r + r )k + cos(a + cr )Nf(r)), (32) 

where k and / are unit orthogonal vectors: A; is a constant vector along the 
string spin and the vector /(r) rotates around it. Introducing constant unit 



vectors i and j, orthogonal to k and to each other we have 

k = J/\J\ = (L 1 + L 2 )/\L 1 + L 2 \, (33) 

i = (L 2 -U)l\L 2 -L x \, j = [k,i\ (34) 

/(r) = A r ~ 1 (— cosa sin(r + r )i + cos(r + r )j), (35) 

N= (l-sin 2 asin 2 (r + r )) 1/2 . (36) 
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The parameters a, ao and tq characterize the initial conditions for the string: 
a is half of the angle between constant vectors of pseudospins L\ and L 2 
having the same length and a and r determine the initial conditions of the 
unit vectors tii, 

Hi = tin sin r — ti i2 cos r, (37) 

in the coordinates i,j, k 

tin = (sin Pi cosa, — cos fy, sin $ sin a), ti i2 = (cos Pi cos a, sin/3j, cos Pi sin a) , 

(38) 

r = (/3i + &)/2, a = (p 1 -p 2 )/2. (39) 

Let us consider the motion of the string in the frame where it is at rest as a 
whole, P — 0. The laboratory time is 

t = x° - x° Q = dr, d = —. (40) 
2ira K J 

— * — * — * 

The string behavior essentially depends on the pseudospins L\ and L 2 , \L\ \ = 
\L 2 \ — L. We have 

m 2 = 8iraL, J = L x + L 2 . (41) 

— * — * 

In case L 1 = —L 2 , the string spin is zero, the string lies in the plane orthog- 
onal to Li and represents a circumference with an oscillating radius, from 
maximal value d to zero and back. 

— * — * 

In the general case, when the angle between L\ and L 2 is in the limits 
< 2a < ii, the value of the string spin is 

TTl 2 

J = cos ex. (42) 

The string represents an ellips with half-axises 

A = dN, B = d sin a cos(r + r ) (43) 

(the large half-axis A is orthogonal to the spin and the small one is parallel 
to the spin) and rotates around the spin with the angular velocity 

\df/dt\ = d- 1 cos aN~ 2 . (44) 

The half-axises are maximal: A — d B = dsina, and the instant angular 
velocity is minimal: c?~ 1 coso;, when the string is in the plane orthogonal to 
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the plane of Lj. Rotating, the ellips shrinks and accelerates. When it reaches 
the plane of Li, it shrinks into a straight-line with half-length dcosa. Its 
angular velocity at this moment is maximal and is equal to the inverse of its 
half-length. Then the ellips expands and slows down and so on. 
In the other extreme case, Li = L 2 , the string spin is maximal 



m 2 



J=ir- ( 45 ) 

(and twice as small as for an open straight-line string of the same mass), the 
string is compressed into a straight-line with half-length d and rotates with 
a constant angular velocity d~ x . 

To quantize the classical solution, one has to know the domain of its 
stability. It is not difficult to show that the above solution is stable for 
all values of the initial conditions except for L = 0. For L = the string 
reduces to a pointlike system with zero mass as it is seen from (27), however 
for arbitrary small L ^ we have system with small mass but with all degrees 
of freedom of the elliptic string. Therefore quantization of our system has no 
meaning for values of pseudospins close to zero (it would lead to a tachion). 



3 Quantization and comparison with experi- 
ment 

Canonical quantization of the model means replacement of the variables 

— * 

P, Z, Li and n« by operators, their Poisson brackets (29) and (30) by com- 
mutators ({, } — > —«[,]) and the constraints following from (27) — by wave 
equations 

(E V / 5-T L ^ 2 -«o)^ = 0, (46) 

(75 - = 0. (47) 

This quantum system is Poincare-invariant because the operators of linear 
and angular momenta, obtained from the corresponding classical expressions, 
satisfy the Poincare algebra [2]. a has been introduced into Eq. (46) to 
account for nonstring interaction at small distances just as it was done for 
quark- ant iquark mesons in Ref. [1]. In general, a could depend on Lj, but it 
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can not increase with L 2 , because the increasing contribution comes from the 
string only (the main assumption of the model). Vacuum fluctuations of the 
higher string modes could also contribute to a . To a first approximation, let 

— * 

us take into account only a part of a which does not depend on L iy that is 
we shall consider constant to be determined from experiment. This 

is the only unknown parameter of the model since the second parameter, the 
string tension a, is known from the analysis of quark- ant iquark mesons [1] 

a = 0,176 ± 0,002 2 . (48) 

Digressing from a factor in the wave function which describes the motion of 
the string as a whole, we can write down the solution of Eqs. (46), (47) in 
the representation where P and Hi are diagonal [2] 

ipjMi(ni,n 2 ) = C{iM]lmi ) lm 2 )Y lmi (n 1 )Y lm2 (n 2 ), (49) 

M=mi+ni2 

where C are the Clebsch - Gordon coefficients. This solution describes the 
internal motion of the quantized string (glueball) with spin j and spin pro- 
jection M and pseudospins I, which take values / = 1, 2, 3, in accord with 
stability of the classical solution considered in the previous Section. The spin 
takes values j = 0, 1, 2, 21. 

We see that the glueballs are space-parity even for all j and I. Since they 
are electrically neutral and charge-parity even, their quantum numbers are 

jGfC = Q+j++ m 

From Eqs. (46) and (47), the glueball mass depends on I only. Introducing 
k = for the leading Regge trajectory, k = 1 for the first daughter trajectory, 
k = 2 for the second one and so on, we can write down for the spin 

j = 21- k. (51) 

Then from Eqs. (46), (47) and (49) the glueball Regge trajectories are equal 
to 

^/(j + k){j + k + 2) = a + —m 2 . (52) 

Glueballs with even spins lie on the leading Regge trajectory, on the second 
daughter trajectory and so on. Glueballs with odd spins lie on the first, third 
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daughter trajectories and so on. All these trajectories j(m 2 ,k) as functions 
of m 2 have at large j the slope 



j'(oo, k) = — |— = 0, 452 ± 0, 005 ~ 2 , (53) 

which is twice as small as for quark-antiquark states and are nonlinear at 
small j. 

The glueball quantum numbers (50) and the mass degeneracy of the states 
with different spins j at the same pseudospin I are remarkable properties of 
glueballs in this model. The lightest states with / = 1 and quantum numbers 
+ ++ , + l ++ and + 2 ++ can be identified with the mesons 

/o(1500), 0+0++, m= 1500 ±10 
/i(1510), 0+1++, m=1518±5 
/ 2 (1565), 0+2++, m = 1542 ± 22 . 

(54) 

These mesons are not the quark-antiquark ones [1] . References to the original 
papers on these mesons and their discussions are given in Ref . [5] , where only 
/o(1500) is considered as firmly established. Let us conservatively estimate 
the lightest glueball mass with I — 1 to be 

mi = 1500 ± 20 . (55) 

This alows one to obtain the constant a from Eq. (52) 

a = 1,81 ±0,04. (56) 

Now, when we know the model parameters we can fix the glueball Regge 
trajectories (52) and to predict masses and quantum numbers of heavier 
glueballs. For / = 2 we have spin-parities and mass 

++ , 1 ++ , 2++, 3 ++ , 4 ++ ; m 2 = 2610 ± 20 . (57) 

For the next glueballs the model predicts 

++ , 1 ++ , 5 ++ , 6 ++ ; m 3 = 3360 ± 25 . (58) 
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From Eq. (52) we get the glueball leading Regge trajectory - - the 
Pomeron trajectory 

j(m 2 , 0) = j(m 2 ) = ^(ao + ^m 2 ) 2 + l - 1, (59) 
with the intersept 

j(0) = 1,07 ±0,03, (60) 

which corresponds to the high-energy data on hadron scattering [6]. Let us 
also note that 

j"(0) = 0,395 ±0,005 " 2 , j"(0) « 0.02 ~ 4 . (61) 

Let us remark in conclusion that glueball decays can be considered within 
this approach if interacting meson fields are introduced instead of meson wave 
functions (second quantization), as it was done in Ref. [7] for open strings. 
Thus, the string quark model provides a single relativistic approach for de- 
scription of quark-antiquark and glueball meson states including prediction 
of the pomeron Regge trajectory, with wide area of critical comparison with 
experiment. 

From viewpoint of this model it is important to make more precise the 
experimental status and properties of the /i (1510) and /2(1565)-mesons and 
to obtain experimental information on possible glueballs in the mass region 
of 2600 and 3300 MeV. 

The author is grateful to V. A. Petrov for stimulating discussions. 
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